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Introduction 

What type of computations can neuronal networks per- 
form? One could reasonable argue that neuronal networks 
should be able to solve "all" computational problems of 
interest, as our brain is made of those networks. However, 
finding specific non-trivial and interesting examples of 
neuronal networks that solve a particular computational 
problem has been quite difficult, remarkably. There are 
very successful examples, though, such a Hopfield net- 
works for memory storage and retrieval, but in particular 
the computational capabilities of spiking networks have 
been far less investigated. Finding new examples is impor- 
tant for a better understanding of the functioning of 
the brain. 

Results 

We find that a family of quadratic programming (QP) 
problems with linear constraints can be solved exactly by 
networks of integrate-and-fire networks, with the only 
approximation being finite number of spikes. We show 
that a network of integrate-and-fire neurons can encode 
an input vector by approximating it as a linear combina- 
tion of stored feature vectors weighted by non-negative 
coefficients. Therefore, the network is able to solve a QP 
problem with non-negativity constraint on the coeffi- 
cients. While in previous rate-based implementation of 
QP problems impose the non-negativity of the firing rate 
artificially in the dynamical system [1], a neuronal net- 
work of interacting neurons satisfies the non-negativity 
constraint for free. We show that a LI and L2 priors of 
the coefficients of the input vectors are encoded in the 
activity of the network as a constant negative current and 
a higher hyperpolarizing reset, respectively. We also 
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show that these networks in the presence of probabilistic 
synapses sample the space of solutions of the QP pro- 
blem, and that this sampling obeys contrast-invariant 
properties. Finally, we find that when the feature vectors 
are dense, the dynamics of the networks have very slow 
modes, which generate slow transients and as a result 
large spiking variability. Despite this large variability, the 
representation of the stimulus is very stable over time. 
Even in the presence of dense features, LI regularization 
reduces spiking variability and allows a very fast conver- 
gence to the best solution. When features are not dense, 
convergence to the best solution is fast regardless of the 
regularization. 

Conclusions 

We have designed a spiking network that can solve 
exactly a QP problem with linear constraints. The net- 
works are able to decompose an input vector into a lin- 
ear combination of stored features with non-negative 
coefficients. This type of networks can be important for 
applications such as odor identification and classifica- 
tion, and biologically-inspired memory storage. 
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